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ABSTRACT 

Two years of microwave background observations with the Cosmic Background Imager (CBI) have 
been combined to give a sensitive, high resolution angular power spectrum over the range 400 < i < 
3500. This power spectrum has been referenced to a more accurate overall calibration derived from the 
Wilkinson Microwave Anisotropy Probe. The data cover 90 deg^ including three pointings targeted for 
deep observations. The uncertainty on the i > 2000 power previously seen with the CBI is reduced. 
Under the assumption that any signal in excess of the primary anisotropy is due to a secondary 
Sunyaev-Zeldovich anisotropy in distant galaxy clusters we use CBI, Arcminute Cosmology Bolometer 
Array Receiver, and Berkeley-IUinois-Maryland Association array data to place a constraint on the 
present-day rms mass fluctuation on 8ft.~^Mpc scales, erg. We present the results of a cosmological 
parameter analysis on the £ < 2000 primary anisotropy data which show significant improvements 
in the parameters as compared to WMAP alone, and we explore the role of the small-scale cosmic 
microwave background data in breaking parameter degeneracies. 

Subject headings: cosmic microwave background — cosmological parameters — cosmology: observa- 
tions 



1. INTRODUCTION 

The Cosmic Background Imager (CBI) is a planar syn- 
thesis array designed to measure cosmic microwave back- 
ground (CMB) fluctuations on arcminute scales at fre- 
quencies between 26 and 36 GHz. The CBI has been op- 

■ crating at its site at an altitude of 5080 m in the Chilean 
Andes sin ce late 1999. Prev io us results have been pre- 

■ sented bv 'Padin et^aP (|200lD . iMason et all l|2003l^ . and 
] ^earson ct al. ( 200J||). The principal observational results 

of these papers were: (i) the first detection of anisotropy 
on the mass scale of galaxy clusters — thereby laying a 
I firm foundation for theories of galaxy formation; (ii) the 
clear delineation of a damping tail in the power spectrum, 
best seen in the mosaic analysis of Pearson ct al.; (iii) 
the first determination of key cosmological parameters 
from the high-£ range, independent of the first acous- 
tic peak; and (iv) the possible detec tion, presented in 
the deep field analysis of lMason et all of power on small 
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angular scales in excess of that expected from primary 
anisotropics. The interpr etation of th ese results has been 
discussed bv lSievers et al. (,20.03) and iBond et al.. i^QM^. 
The CBI data, by virtue of their high angular resolu- 
tion, were able to place constraints on cosmological pa- 
rameters which are largely independent of those derived 
from larger-scale experiments; for instance, 10% mea- 
surements of iltot and Us using only CBI, DMR and a 
weak Hq prior. The small-scale data also play an impor- 
tant role in improving results on certain key parameters 
(ilfc/i^, ns, Tc) which are less well-constrained by large- 
scale data. 

Theoretical models predict the angular power spec- 
trum of the CMB 

Ci = (|a,„p) (1) 
where the agm are coefficients in a spherical harmonic 
expansion of temperature fluctuations in the CMB, 
AT/TcMB, where Tcmb « 2.725 K is the mean tem- 
perature of the CMB, and the angle brackets denote 
an ensemble average. These theories also predict a se- 
ries of acoustic peaks in the angular power spectrum 
on scales < 1° > 200), and a decline in power to- 
wards higher £ due to photon viscosity and the thickness 
of the last scattering surface. Earl y indications of th e 
flrst acoustic peak were presented bv lMiller et"a D ('1999f); 
definitive measure ments o f the first a nd se c ond peaks 
were r eported bv Ide Be rnardis ct al. (2 00(11). iLee et alJ 
(l200l.lNetterfield et alt (^2002.') . .Halverson et alJ (|2002|) . 
IScott et alJ (120031) . and IGrainge et all (120031) ^^. The 
last of these experiments reached £ ^ 1400. The CBI 
(Padin et al. 2002) has complemented these experiments 
by covering an overlapping range of £ extending to 
£ ^ 3500. The Arcmin ute Cosmology B olometer Ar- 
ray Receiver (ACBAR) (|Kuo et al.l 120041) has recently 

In the parameter analysis of §4 we use the latest VSA data 
IDickinson et alJ2004l) . which was released shortly after this paper 
was first submitted. 
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covered a similar range of I as the CBI at higher fre- 
quency; the Berkeley-Ilhnois-Maryland Association ar- 
ray (BIMA) has also made 30 GHz measurements at 
£ ~ 5000 which probe the secondary Sunyaev- Zeldovich 
effect (SZE) anisotropy f Pawson et'a^l2002^ . These 
experiments — which employ a wide variety of instrumen- 
tal and experimental techniques — present a strikingly 
consistent pic ture whi ch supports inflationary expecta- 
tions (sec Bon d etal] 12002 for a review). However the 
results at intermediate angular scales (500 < I < 2000) 
currently have comparatively poor £-space resolution, 
and the high-£ results are difficult to compare conclu- 
sively owing to the low signal-to-noise ratio (~ 2-4). The 
results presented here improve the situation by: (i) ex- 
panding the coverage of the CBI mosaics for higher £ 
resolution, (ii) integrating further on the deep fields, and 
(iii) combining the deep and mosaic data for a single 
power spectrum estimate over the full range of i covered 
by the CBI. 

The CBI resuhs presented bv lMason et all f2003) and 
IPears on et al- (2003-) were based on d ata obta ined be- 
tween January and December of 2000. Mas on et alJ an- 
alyzed the data resulting from extensive integration on 
three chosen "deep f ields" to const rain the small-scale 
signal; the analysis of ! Pearson et alJ used data with shal- 
lower coverage of a larger area ( "mosaics" ) to obtain bet- 
ter Fourier-space resolution. Further observations were 
conducted during 2001; these were used to extend the sky 
coverage of the mosaics in order to attain higher resolu- 
tion in £, and to go somewhat deeper on the existing deep 
fields. This paper presents the power spectrum resulting 
from the combination of the full CBI primary anisotropy 
dataset, which comprises data from years 2000 and 2001 
on both mosaic and deep fields. Two of the mosaic fields 
(14 h and 20 h) include deep pointings; there is also a 
third deep pointing (08 h), and a third mosaic (02 h). 
The CBI data have been recalibrated to a more accu- 
rate power scale derived from the Wilkinson Microwave 
Anisotropy Probe [WMAP). 

The organization of this paper is as follows. In §[51 we 
discuss the observations and l/FAL4P-derived recalibra- 
tion. In §121 we present images and power spectra derived 
from the data and explain the methodology employed in 
their derivation. In §^we use these results to constrain 
cosmological parameters based on standard models for 
primary and secondary CMB anisotropics. We present 
our conclusions in §0 

2. OBSERVATIONS AND CALIBRATION 

The analysis in this paper includes data collected in 
the year 2001 in addition to the year 2000 data previ- 
ously analyzed. In January through late March of 2001 
there was an unusually severe "Bolivian winter" which 
prevented the collection of useful data; regular observa- 
tions resumed on 2001 March 28 and continued until 2001 
November 22. The weather in the austral winter of 2001 
was considerably less severe than it had been in 2000, 
so that significantly less observing time was lost due to 
poor weather. 

In 2001 we concentrated primarily on extending the 
mosaic coverage in three fields in order to obtain higher 
resolution in t. We also made a sma ll number of observa - 
tions in the deep fields discussed bv lMason et alJ ll2003f). 
The original field selection is discussed bv iMason^taLl 



lj2003(l and IPearson et alJ l|2003|) . Since our switching 
strategy employs pairs of fields separated in the east- 
west direction, contiguous extensions were easiest in the 
north-south direction. The extensions to these fields were 
selected to minimize point source contamination. In two 
cases (the 02 h and 14 h fields) this procedure resulted 
in extensions to the north, and in one case (20 h) an 
extension to the south. The images for the combined 
2000-f2001 mosaic observations are shown in Figure ^ 
and the sensitivity maps are shown in Figure (21 The to- 
tal areas covered are 32.5, 3.5, 26.2, and 27.1 deg^ for the 
02 h, 08 h, 14 h, and 20 h fields respectively". 

The CBI two-year data were calibrated in the same 
manner as the first-year data |Mason ct al. 2003), ex- 
cept that the overall calibration scale has been adjusted 
in light of the r ecent WMAP observations of Jupiter 
ijPaee et al.ll2003() . The flux density scale for the first- 
year data was determ ined from th e absolute calibra- 
tion measurements by iMason et al.l (|1999fl which gave 
a temperature for Jupiter at 32 GHz of Tj = 152 ± 5K 
(note that all planetary temperatures discussed in this 
paper are the Rayleigh- Jeans brightness temperature of 
the planet minus the Rayleigh- Jeans temperature of the 
CMB at the same frequency) . Page ct al. ( 2003) have de- 
termined Tj (32 GHz) — 1 46. 6 ± 2. K from measurements 
relative to the CMB dipole. Figure 01 sh ows measure- 
ment s of Jupiter with the CBI on the old ({Mason et alJ 
Il999|) calibration scale, as well as the WMAP measure- 
ments. The slopes of the CBI and WMAP spectra are 
in agreement to better than Icr, and the 32 GHz val- 
ues are also within Ict; these results support the original 
CBI amplitude and spectral slope calibrations. Since the 
WMAP and Mason et al. measurements are independent 
we adopt a weighted mean of the two and base the CBI 
calibration on rj(32GHz) = 147.3 ± 1.8 K. This 3% 
reduction in the CBI flux density scale corresponds to 
an overall scaling down of the CBI power spectrum by 
6% in power. This scaling is consistent with the original 
3.5% flux density scale uncertainty (7% in power). The 
new CBI calibration has an uncertainty of 1.3% in flux 
density (2.6% in power). 

3. DATA ANALYSIS 

The basic methods of CBI data a nalysis and sp ectrum 
extraction are described fully by IMason et al.l (12003), 
IPearson et all l|2003t) . and lMvers et all l|2003'). The pri- 
mary differences in this analysis are: an improved esti- 
mate of the thermal noise, which has allowed us to bring 
the mosaic data to bear on the "high-£ excess" evident in 
the deep data; new data cuts needed to deal with point 
sources in the mosaic extensions; and a revised ^-binning 
appropriate to the expanded sky coverage and variable 
noise level of the new data. These aspects of the anal- 
ysis, and the resulting power spectrum, are described in 
the following subsections. 

3.1. Thermal Noise Estimates 

We estimate the thermal noise variance for each uv 
data point in each scan of the CBI dataset by the mean 
squared-deviation about the mean, and subsequently use 
a weighted average to combine the estimates from differ- 
ent scans. It is necessary to make a small correction to 

These are the areas, counting LEAD and TRAIL fields sepa- 
rately, mapped to an rms sensitivity of lOmJy/beam or better. 
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Fig. 3. — Comparison of Jupiter temperatures measured with the 
CBI and with WMAP. The spectrum of Jupiter in this frequency 
range is not thermal owing to an absorption feature. The individual 
channel CBI temperatures of Jupiter in the 26-36 GHz range are 
shown by the filled blue circles, and those of WMAP in the range 
22—94 GHz are shown by the open pink squares. The dotted and 
solid lines, respectively, show the best fit slopes to WMAP and 
CBI data over the 20 to 40 GHz range. The CBI temperatures 
shown here were d etermined assuming Tj(32GHz) = 152 ± 5K 
JMason et al.ll999l) : the systematic error bar on this calibration is 
shown as a dotted red line. The WMAP error bars include both 
random and systematic errors. 



the estimated variance for each {u, v) data point due to 
the finite n umber of samples which enter the estimate. 
iMason eT^ l. (2003) present simple analytic arguments 
placing this correction at 8% in variance, and estimate 
a 2% uncertainty in the variance. We have since im- 
proved our estimate of the CBI thermal noise variance 
resulting in a variance correction factor of 1.05 ± 0.01. 
This is 1.5(7 from the factor (1.08 ± 0.02) applied to the 
year 2000 data; the overestimate of noise in Mason et al. 
will have caused a slight underestimate — by 42 /iK^ — of 
the ex cess power level at high £ in the original CBI deep 
result ijMason et al.ll2003|) . 

We calculated the noise variance correction in several 
ways: a first-order analytic calculation of the noise dis- 
tribution; a numerical (FFT-based) integration of the 
distribution; and analysis of simulated data. The simula- 
tions use the actual CBI data as a starting point, replac- 
ing each 8.4 second integration with a point of zero mean 
and drawn from a Gaussian distribution with a dispersion 
derived from the estimated statistical weight of the data 
point. This method accounts for variations in the statis- 
tical weight from baseline to baseline and variations in 
the number of data points per scan due, for instance, to 
rejected data. This gives a result of 1.044 ±0.002 (statis- 
tical error). The FFT calculation agrees within the 0.2% 
uncertainty in the variance. The first-order analytic cal- 
culation is lower by 0.6% but should be considered to be 
only a crude chec k. Deta i ls of t he noise variance estimate 
are presented bv lSiever"^ 1)20031) . 

The simulations were analyzed entirely in the visibil- 
ity domain. As a further check on the variance cor- 
rection calculations, we gridded the visibility data fol- 
lowin g the procedure u sed in analyzing the real CBI 
data llMvers et alJl2063|) . Monte Carlo calculations of 
the of the gridded estimators (using the inverse of 
the full noise matrix) yield a noise correction factor of 
1.050 ± 0.004 (statistical error). This is why we adopt a 
value of 1.05, since it is the gridded estimators that are 



used in the power spectrum estimation. We conserva- 
tively adopt a 1% uncertainty in the variance correction 
although the level of agreement between different meth- 
ods of estimating this is ~ 0.5%. 

The importance of the improved accuracy of the ther- 
mal noise calculation is illustrated by considering the 
year 2000 CBI data a.t £ > 2000. Referenced to the cur- 
rent calibration and noise s cale, t he ye ar 2000 mosaic 
data presented by Pearson et al.l l|2003D yield a band- 
power of 206 ±178 ^K^. We considered combinin g these 
data with the year 2000 CBI deep field data ( Maso n~t alJ 
|2003|) . The thermal noise variance in this last bin for the 
year 2000 mosaic, however, is 4307 /zK^, which yields a 
86 /iK^ systematic uncertainty due to the thermal noise 
variance correction alone given the previous 2% uncer- 
tainty in the noise variance. The thermal noise variance 
in the year 2000 deep field data is 1293 /iK^, resulting 
in a 26 /iK^ uncertainty due to the noise estimate; this 
was substantially less than the greatest systematic un- 
certainty, the residual point source correction at 50 /iK^. 
It was clear that the mosaic data would contribute lit- 
tle to our understanding of the signal at ^ > 2000. In 
contrast, the thermal noise variance for the 2000+2001 
mosaic±deep data in the last {£ > 1960) bin is 2054 /iK^, 
which with our present 1% knowledge of the noise vari- 
ance results in a contribution to the systematic error bud- 
get slightly lower than that from noise in the year 2000 
deep data, and subdominant to the residual point source 
contribution. 

3.2. Treatment of the Discrete Source Foreground 

The treatment of the discrete source foreground is 
si milar to that adopt ed in the earlier CBI analyses 
of IMason etldl (|2003(1 and IPearson et a l. (2003). All 
sources above 3.4 mJy i n the 1.4 GH z NRAO VLA 
Sky Survey (NYSS: iCondon et all 1199^ were included 
in a constraint matrix and projected out of th e data 
lIBond. .Taffe. fc KnoxllT9M IHalverson eTaLllMil . This 
is roughly equivalent to completely downweighting the 
synthesized beam corresponding to each of these sources 
and effectively eliminates 25% of our data. We correct 
for sources below the 3.4 mJy cutoff in NVSS statistically. 
The statistical correction reduces the power in the high-i? 
bin by ^ 20% (see Figure 0J. We have also obtained in- 
dependent 30 GHz measurements of the bright sources in 
the CBI fields with the OVRO 40-m telescope. For pre- 
sentation purposes we subtract these flux densities from 
the maps with reasonable results, although some resid- 
uals are visible. The constraint matrices eliminate the 
impact of any errors in the source subtraction, and the 
power spectrum results are unchanged even if no OVRO 
subtraction is performed. 

Although the extensions were chosen to minimize point 
source contamination, the larger size of the expanded 
mosaics and the requirement that the extensions be con- 
tiguous with the already highly-optimized original CBI 
mosaic fields resulted in a handful of sources brighter 
than those present i n the year 2000 CBI mosaic data 
ijPearson et al.l l2003|) . In particular the 02 h extension 
contains the Seyfert galaxy NGC 1068 (S'soghz ~ 0.4 Jy) 
which we found was not effectively removed by the con- 
straint matrix. To deal with this we excluded CBI point- 
ings around this source (as well as one other bright source 
in the 02 h field, and one in each of the 14 h field and 20 h 
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fields) until the maximum signal-to-noise ratio on any 
discrete source in the total mosaic areas — before subtrac- 
tion or projection — was less than some threshold X. In 
our final analysis we have adopted X = 50, eliminating 9 
out of 263 CBI pointings. To check this we analyzed the 
data with a more stringent SNR cutoff of X = 25 and 
found no significant change in the power spectrum. 

3.3. Power Spectrum Analysis 

The dataset combines very deep pointings (and thus 
low noise levels) on a few small areas with substantially 
shallower coverage of wider areas. The signal at low- 
(. is stronger and the features in the power spectrum 
are expected to be more distinct, so we seek to use the 
wider coverage for maximum i resolution on large angu- 
lar scales. Most of the statistical weight in the dataset at 
small scales comes from the deep integrations, and since 
the sky coverage of these is quite limited the ^ resolution 
is lower. In this regime the signal-to-noise ratio is lower 
and we seek to compensate by combining many Fourier 
modes. In order to present a single unified power spec- 
trum which makes use of information from all the data 
over the full range of angular scales we adopt bins which 
are narrowest at low £ {51 = 100), and increase in steps 
towards a single high I bin above I ~ 2000. The bin 
widths were chosen to yield maximum t resolution while 
keeping typical bin-to-bin anti-correlations to less than 
~ 30%. We have chosen two distinct binnings of the data 
which we call the "even" and "odd" binnings. These bin- 
nings are not independent but are helpful in determin- 
ing whether particular features are artifacts of the bin 
choice. Subsequent statistical analyses — including pri- 
mary anisotropy parameter estimation and the cts anal- 
ysis of secondary anisotropy — employ only one binning 
(the "even" binning). 

The updated CBI power spectrum is shown in Fig- 
ure 21 and tabulated in Table ^ Results are presented 
in terms of bandpower (AT^ = T'^MB^i^ + l)C£/27r), 
which is assumed to be flat within each bin; also shown 
are the values of the noise power spectrum Xb- This 
table presents both "even" and "odd" binnings of the 
CBI data. Window functions for the two binnings are 
presented in Figures El and The procedures for cal- 
culating w indow functions an d noise power spectra are 
defined bv lMvers et"al] l|2003() . 

The possible detection of power in excess of 
the expected pr imary anisotropy at high-€ by 
Mason et al.l 112003) has s t irred considerable interest (e.g. 
Koniatsu fc SeliaH l2002t lOh. Coorav. fc Kamionkowskil 
20031 i Subramanian. Sesha dri. fc Barrow! 120031: 
Griffiths. K unz. & Silk 2003), and the results we 
present in this paper improve the bandpower constraint. 
Binning aU the 2000-1-2001 data above I = 1960 together, 
we find a bandpower of 355 ± 103 /zK^ (random error 
only). By way of comparison, the IMason et alJ result, 
referenced to the current CBI calibration scale and noise 
correction, is 511 ± 156 /iK^; the new result is thus 
25% lower but within ~ la of the Mason et al. result, 
although the datasets are not independent. Table (21 
presents the (. > 2000 bandpower constraints from these 
and three other combinations of the full CBI dataset, 
all referenced to the current calibration scale and with 
our best noise variance estimates. For purposes of 
comparison this table shows only the random errors 



derived from the Fisher matrix at the best fit point 
(which includes couplings to other bins); in addition 
there is a common overall uncertainty of 48 /^K^ from 
the discrete source correction. 

In order to constrain the excess more accurately we 
have explicitly mapped the likelihood in the last bin, al- 
lowing for the following systematic errors in the anal- 
ysis: uncertainty in the statistical source correction 
(48 /iK^); uncertainty in the thermal noise power spec- 
trum (20 /iK^); and the 56 /zK^ dispersion in the high-£ 
bandpower caused by the uncertainty in the bandpower 
of the neighboring bin^"*. We determine confidence in- 
tervals on AT^ of 233-492^K2 (68%) and 110-630AtK2 
(95%). From the 68% confidence hmit we can state our 
result as 3551^22 mK^. This result is consistent with 
but lower than that derived from the earlier analysis of 
CBI deep fields; and while the detection of power re- 
mains statistically significant, the detection of power in 
excess of the band-averaged power expected from the pri- 
mary anisotropy 80-90 ^K^) is marginal. A slightly 
more significant detection is obtained by combining CBI, 
ACBAR, and BIMA data, and we present the results of 
such an analysis in S 14.31 

We have also computed the value of the high i bin 
for several statistically independent splits of the total 
(2000-1-2001 deep plus mosaic) dataset. In all cases the 
power spectra are consistent. The most sensitive of these 
splits was a division of the dataset into two halves by 
field (02 h plus 08 h, and 14 h plus 20 h), in which case 
the high £ bandpowers were within 1.3(t of each other 
and consistent with the best value of 355 ^K^. 

4. INTERPRETATION 

4.1. Basic C'osmological Parameters from the Primary 
Anisotropy 

We use a modified version of the publicly available 
Markov Chain Monte Carlo (MCMC) package COS- 
MOMC^^ ijLewis fc Bridle 2002) to obtain cosmological 
parameter fits to the CMB data. The code has been 
tested extensively against our fixed resolution grid based 
method, which w e applied to t he firs t yea r CBI data 
in the p apers bv iSievers et alJ f2003D and iBond et alJ 
l|2004D . IBond. Contaldi fc PogosvanI l|2003t) show that 
the agreement between the two methods is good when 
identical data and priors are used. Advantages of the 
MCMC method include a reduced number of model 
spectrum computations required to accurately sample 
the multi-dimensional likelihood surfaces and automatic 
rather than built-in adaptivity of the parameter sets sam- 
pled. 

Our typical run involves the calculation of 8 Markov 
chains over the following basic set of cosmological pa- 
rameters: LOh = ^bh^, the physical density of baryons; 
LOc = ^ch"^ , the physical density of cold dark matter; f2A, 
the energy density in the form of a cosmological con- 
stant; rtg, the spectral index of the scalar perturbations; 
As, an amplitude parameter for the scalar perturbations; 
and Tc, the optical depth to the surface of last scattering. 

The error quoted from the Fisher matrix at the best fit point, 
±103 /iK^, includes this contribution, and it is only added in sepa- 
rately here because the likelihood mapping procedure keeps other 
bins fixed. 

|http : //cosmologist ■ inf o/cosmomc/ 1 
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Fig. 4. — The 2000+2001 CBI Spectrum. The "even" binning is shown in red and the "odd" binning in hght blue. Orange 
stars indicate the thermal noise variance; green triangles indicate the statistical source correction which has been subtracted from the 
power spectrum. The solid black line is the WMAP ACDM model with a pure power-law primordial spectrum (model spectrum is file 
wmap_lcdjn_pl_model_yrl_vl.txt, available on the WMAP website ^http:/ /lambda. gsfc.nasa.gov J . 



TABLE 1 
CBI Bandpower Measurements 





Even Binning 




Odd Binning 




Bin 


^-Hange 


Bandpower 




£-Hange 


Bandpower 














(mK2) 


(mK2) 


1 


0-300 


7091 ± 1882 


3176 


0-250 


7860 ± 4151 


8196 


2 


300-400 


2059 ± 717 


489 


250-350 


4727 ± 1157 


796 


3 


400-500 


1688 ± 457 


377 


350-450 


961 ± 454 


397 


4 


500-600 


2415 ± 545 


449 


450-550 


2369 ± 504 


390 


5 


600-700 


1562 ± 391 


423 


550-650 


2081 ± 480 


455 


6 


700-800 


2201 ± 490 


577 


650-750 


1494 ± 400 


460 


7 


800-900 


2056 ± 436 


631 


750-850 


2346 ± 476 


582 


8 


900-1000 


1158 ± 396 


743 


850-950 


2117 + 482 


770 


9 


1000-1140 


797 ± 275 


674 


950-1070 


305 ± 239 


636 


10 


1140-1280 


780 ± 263 


726 


1070-1210 


1226 ± 300 


694 


11 


1280-1420 


586 ± 278 


933 


1210-1350 


423 ± 269 


820 


12 


1420-1560 


1166 ± 361 


1064 


1350-1490 


1020 ± 333 


1040 


13 


1560-1760 


196 ± 223 


941 


1490-1660 


714 ± 279 


960 


14 


1760-1960 


-4 ± 203 


386 


1660-1860 


-98 + 201 


834 


15 


1960+ 


355 ± 103 


2184 


1860-2060 


243 + 229 


457 


16 








2060+ 


346 + 113 


2385 



Each chain is run on a separate 2-CPU node of the CITA 
McKenzie Beowulf cluster for a typical run-time of ap- 
proximately 9 hours when the proposal densities are es- 
timated using a previously computed covariance matrix 
for the same set of parameters. The chains are run un- 
til the largest eigenvalue returned by the Gelman-Rubin 
convergence test reaches 0.1. We run the chains at a tem- 
perature setting of 1.2 in order to sample more densely 
the tails of the distributions; the samples are adjusted 



for this when analyzing the chains. 

All of our parameter analysis imposes a "weak- 
ft," prior comprising limits on the Hubble parameter 
(45kms"iMpc"i < Hq < 90 kms^^Mpc^^) and the 
age of the universe (to > lOGyr). We primarily consider 
flat models (f2tot = 1) in this work, and unless other- 
wise stated a flat prior has been imposed. Within the 
context of flat models the weak-h prior influences the re- 
sults very little. We include all of the bandpowers shown 
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Fig. 5.— The 2000+2001 CBI window functions ("even" binning). 
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Fig. 6.— The 2000+2001 CBI window functions ("odd" binning). 



TABLE 2 
Comparison of high-^ results for 
different subsets of cbi data 



Dataset 


Bandpower (/jK^) 


2000+2001 Deep+Mosaic 


355 ± 103 


2000+2001 Deep 


393 ± 134 


2000 Deep 


514 ± 158 


2000 Deep+Mosaic 


457 ± 122 


2000 Mosaic 


206 ± 178 



Note. — Results for the high-£ bin on the 
WMAP power scale, with current noise correc- 
tion applied. For further discussion see text. 



in Tabled except for the highest and lowest I band. The 
highest band is excluded due to possible contamination 
by secondary anisotropics; the first band is poorly con- 
strained and provides no useful information. 



In Table 01 we compare the constraints obtained when 
including only the WMAP measurements with those 
obtained when also including the new CBI bandpow- 
ers and a compilation of "ALL" present CMB data^^ 
for the weak-/i prior case. We include both total in- 
tensity and TE spectra from WMAP in our analysis. 
For Boomerang and ACBAR, recalibrations and their 
uncertainties were applied using t he power spectrum 
based method described in Bon d. C'ontaldi fc PogosvarJ 
|)200.1D which obtains maximum likelihood calibration 
parameters as a by-product of the optimal, combined 
power spectrum calculation with multiple experiment re- 
sults. De tailed results for the fits arc summarized in 
Table 1 of iBond. Contaldi fc PogosvanI (|2D03). We note 
that this method gives calibrations in agreement with 



1^ WMAP IBennett et al."'2003'); VSA f Dickinson et al."2004D; 
DASI iHalverson et al. 2002); ACBAR (Kuo et al. 2004); MAX- 
IMA <Lee et abi2001.) : and BOOMERANG iRuhl et al 200a) . 
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TABLE 3 

COSMOLOGICAL CONSTRAINTS FROM THE "WMAP ONLY", "CBI 

+ WMAF\ AND "CBI + All" for an assumed f2tot = 1.0. 



Parameter 


WMAP only 


CBI + WMAP 


CBI + ALL 




U.UZ'iO_Q QQJg 


0225+0-0011 


n 0995+0 0009 




Q ,23+0 01'' 


107+0-012 
^■l'-"-0.013 


Q -,+0.010 
'-'■^^^-0.009 


Ua 


7, +0.08 
^■'^-0.08 


n 7c:+0.05 
^- '"-0.05 


74+0.05 


TC 


n 10+0. 03 

1 01+005 


13+0-02 


11+0-02 
'^-11-0.03 


ris 


96+0-03 
^-^"-0.03 


95+0-02 


lOio^s 


27 7+5.5 


90 9+2.8 
. ^ 2 


21.9+23 


Ho 


72 1 +^-* 
'^•-'--5.8 


73.4lt'? 


71.911-1 


Age (Gyr) 


1 o q+0.3 
J^J-J_0.3 


-TO 7+0.2 
l'J-'-0.3 


io 7+0.2 




0.29to.08 


^-^'-'-0.05 


26+0-04 


0-8 




n 70+0. 08 


P on+0.06 



Note. — We included weak external priors on the Hubble 
parameter{45 km s~i Mpc~i < Hq < 90 kms~lMpc~l) and 
the age of the universe (to > 10 Gyr). The flatness prior has the 
strongest effect on the parameters by breaking the geometrical 
degeneracy and allowing us to derive tight constraints on Hq 
and Hm- 



TABLE 4 

COSMOLOGICAL CONSTRAINTS FROM "CBI + All" 
DATA PLUS LSS CONSTRAINTS 



CBI + ALL + 2df CBI + ALL + LSS 



^6/12 


0224+0-0008 


n 0995+0-0009 




Q ,,7+0.007 

'-'■^^'-0.006 


n ,,0+0.007 


Ua 


7, +0.03 
'1-0.03 


Q 7, +0.03 
^-'1-0.04 


TC 


10+0-02 
'-'■ll'-0.02 


11+0-02 
'-'-11-0.0 


ris 


95+0-02 

'-'-^'-'-0.02 


95+0-02 

'-'-^'-'-0.02 


lOiOAs 


21 6+2-1 


22.31^-1 


Ho 


69.6+2-5 


69.6+29 


Age (Gyr) 


1 o 7+0.2 


, 7+0.2 




29+0-03 
'-'-^^-0.03 


29+0-0-1 
^-■^^-0.03 


0-8 


09+0.05 


r, oq+0.05 


Note. — 


The priors are the 


same as in Table 01 In 



addition we have added a LSS prior in the form either of 
constraints on the combination a^Vl^^ and the shape 
parameter Fcfl, or using the 2dfGRS power spectrum 
results. 



TABLE 5 

COSMOLOGICAL CONSTRAINTS INCLUDING A RUNNING SPECTRAL INDEX 





WMAP + LSS 


CBI + WMAP + LSS 


CBI + ALL + LSS 




0249+0-0025 
u.uz4y_Q 0025 


n 0999+0-0019 

U.UZZZ_g 00^7 


0918+0-0013 
O.U2i8_o 0014 




116+0-013 
U.11D_Q 0J3 


120+0-013 
U.iZU_0 0j^2 


124+0-011 
U.iZ4_Q 011 


Ua 


74+0-08 
'^-0.08 


69+0-08 


67+0-0'' 
'-'■O'-0.07 


■TC 


r, 09+0.08 
^^-■^^-0.07 


24+0-05 
1-^^-0.07 


21+0-03 


ris 


-1 n+0.08 
l-l'-0.08 


90+0-06 


0.8810-05 


as 


-0 061+0-037 
U.UDi_Q 037 


-0 085+0-031 

U.UOO_0 030 


-0.08710-028 


lOiOyls 


36 1+10-0 
■^''•1-10.0 




28.11^:3 


Ho 


75.71^-27 


68.9lg-2 


67.0+5-1 


Age (Gyr) 


iq 9 + 0.5 
l''-^-0.5 


10 7+0:3 
l'^- '-0.4 


13.8lo;2 




26+0-08 
U.ZD_Q 0g 


r. n,+0.08 
""Jl-O.OS 


Q OO+0.07 

'^■■'•-'-0.07 


0-8 


1 0+0-1 


92+0-08 
'J-y^-0.08 


91+0-07 
'-'■yi-0.07 



Note. — Cosmological Constraints including a running spectral index 
as = dus/dlnk obtained from the CMB and our conservative LSS prior. We 
find all combinations prefer a negative value for as with significances above 
the 2(7 level for the combinations CBI + WMAP and CBI + ALL. 
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those obtained from the WMAP/CBl cross-cahbration 
via Jupiter, and a map-based comparison of Boomerang 
and WMAP gives a very similar recahbration and error 
for Boomerang to those used here (E. Hivon 2003, pri- 
vate communication) . The original reported calibrations 
of DASI, Maxima, and VSA were used. Although the 
optimal spectrum procedure also produces best fit values 
with errors for the beam parameters of each experiment, 
we have used the reported beams and their uncertainties 
in each case for the parameter estimates given in this 
paper. 

The "CBI + ALL" parameters and their errors in Ta- 
ble 13 can be compar ed with the "M arch 2003" values 
given in Table 5 of B ond. Contaldi fc Poeosvan (2003). 
These were evaluated using the MCMC method with the 
calibrations for CBI used here, but no recahbration with 
decreased errors for Boomerang and ACBAR. The results 
are quite similar. 

Our main results for the flat plus weak-h case are 
summarized in Figure \7\ showing marginalized one- 
dimensional distributions for the basic six parameters to- 
gether with three other derived parameters: the Hubble 
parameter Hq in units of km s~^Mpc~^, the total age 
to of the universe in Gyr, and the total energy density 
of matter Qm in units of the critical energy density. We 
show three curves for each parameter corresponding to 
the "WMAP only", "CBI + WMAP\ and "CBI + ALL" 
cases. They show how the inclusion of high-£ bandpowers 
is crucial to excluding significant tails in the distributions 
that remain because of the limited £-range of the WMAP 
results. 

Of particular significance is the effect of including the 
CBI band powers on the correlated trio Ug, tq, and w?,. 
The " WMAP only" case shows long tails towards high 
values for the three parameters which are only excluded 
when the CBI or the "CBI -I- ALL" combinations are 
included. We do n ot include a cutoff o n the value of 
Tc as was done bv lSpergel et all l)2003j) . Their cutoff 
has a rather strong effect also on the tails of the distri- 
bution in Us and Wf,. We rely only on the addition of 
extra data. This can be seen in Figure |H1 which shows 
the marginalized distribution in the n^-Wf, plane for the 
" WMAP only" , "CBI + WMAP' , and "CBI + ALL" 
cases. 

The results of the CMB+LSS parameter analyses are 
presented in Table 0] We consider two cases to illus- 
trate the impact of large scale structure observations 
on the cosmological parameter distributions: (i) the 
Two Degree Ga,laxy Redshift Survey (2dfGRS) results of 
iPercival et al.l l)2003D . and (ii) a more conservative LSS 
prior that straddles most weak lensing and cluster re- 
sults for the amplitude ag l)Bond. Contaldi fc Pogosvai] 
120031 and references th erein), but a weaker version of 
the 2dFGRS and SDSS llTegmark et al.ll200l results for 
the shape of the matter power spectrum than the di- 
rect application of the 2dfGRS data gives. Explicitly 
the prior on the amplitude is as^^^^ = 0-47lo'o2 ^o.o8' 
distributed as a Gaussian (first error) convolved with a 
uniform (top-hat) distribution (second error), both in 
\n{as^m )■ The prior on the effective shape parame- 
ter is Feff = 0.21+o °3;+o °g. Again the small-scale CMB 
results substantially improve the constraints in compari- 
son to what is obtained with only the larger angular scale 



CMB data. Figure |51 shows the rc-cg plane, illustrating 
the exclusion of the high values along the line of near- 
degeneracy which results when CBI and ACBAR data 
are added to WMAP+hSS. 

All of the above analysis assumes iltot = 1- It is well 
known that revoking this assumption yields substantially 
worse parameter estimates when CMB data are analyzed 
in isolation (e.g.,'Efstathiou & Bondl'lOOO^'Soergel^et^^D 
|22g3; Bond. Contaldi & Pogosyan 2003; Tegmar kTtaLl 
120041. and references therein). The main parameters af- 
fected are $7^, ^A, and Hq; typically low Hubble pa- 
rameters and larger ages to are favored in this case. For 
CBI-I-ALL we find a factor of ~ 2 — 3 degradation in 
the precision of the constraints on fl\ and Qm- The 
best value for the curvature in this scenario is flk = 
-0.052to;Q3g. Results on fifc/i^, Qch'^, tc, and are 
not significantly affected^^. Thus CMB data alone yield 
a robust determination of the non-baryonic dark matter 
density, and a determination of the total baryon content 
of the universe consistent with th ose derived from deu - 
terium absorption measurements ijKirkman et alJl2003(l . 
as well as limits from other light-element methods (e.g. 
iBania. Rood, fc Balseill2002l and references therein). 

4.2. Running of the scalar spectral index 

Inflation models rarely give pure power laws, with 
ns{k) constant, even over the limited ranges of wavenum- 
ber k that the CMB-I-LSS data probe. In most models, 
the breaking is rather gentle, with small corrections hav- 
ing been entertained since the early eighties. Much more 
radical forms for ns{k) are possible. The gentle form 
most often adopted involves a running index described 
by a logarithmic correction: 

dlnP T f k \ 

ns[k) = -j^—j- = ns(koj + asln [— j , (2) 

where as — dus/dhik. Here P{k) is the primordial post- 
inflation power spectrum for scalar curvature perturba- 
tions and ko is a pivoting scale about which ns{k) is 
expanded. The effect is that for negative as the slope is 
flattened below ko and steepened above ko, i.e., power is 
suppressed on scales both greater than and less than feg. 

There has been much focus recently on whether 
the data require such a running index, motivated 
by the incorporation of Lym an-g absorpti on data 
in the WMAP analyses of iSpergel et all_ (|2003), 
and explored further by, e.g. , Bridle et al. (2003), 
[Bas tero-Gil. Fre ese. fc M ersini-Houghton (2003), and 
iMukhcriec fc Wangl (20031. Bond. Contaldi fc Po gosvanI 
l|2003i) have shown that the CMB data marginally favor 
a non-vanishing negative running term. The effect is 
driven by the requirement to reconcile an apparent lack 
of power on the largest scales observed by WMAP with 
observations on arcminute scales such as those reported 
in this work. In this regard, CBI adds a significant lever 
arm beyond WMAP to higher £, and the CBI/ WMAP 
cross-calibration presented here therefore helps to fur- 
ther constrain the allowed variation of ns{k). 

Table [SI shows the parameters we obtain when our ba- 
sic parameter set is expanded to include a running term 

This explains the mechanism for degraded estimates of other 
parameters: increased uncertainty in Hq coupled with fixed values 
of flfjh'^ and Qch^ leads to the increased uncertainty in f2m, also 
causing an increased uncertainty in Q/^. 
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Fig. 7. — Marginalized likelihood curves for a range of individual cosmological parameters, each shown for three CMB datasets: " WMAP 
only" (blue/dotted); "CBI + WMAF^ (red/dashed); and "CBI + ALL" (green/solid). In all cases a flat plus weak-/i prior is used. 



as = drig/dliik, with the LSS prior apphed for the three 
cases. We have not hmited as by any theoretical prior 
prejudices, but have just allowed it to vary over the range 
—0.2 < as < 0.2. The final 1-d marginalized distribu- 
tions for a number of combinations of data and priors 
are shown in Figure [TUl Analyzing the WMAP data 
alone, we find as = — 0.077lg Qgg; including the CBI re- 
sults favors a more negative value ag = — O.lOSlQggg. 
Adding LSS constraints reduces the uncertainties some- 
what, yielding Ug — — O.OSS^Q QgJ. Estimates for the 
optical depth tc and linear amplitude erg are generally 
higher and those for the spectral index at the chosen 
pivot scale ns{ko = 0.05 Mpc~^) are lowered. Figure ITTI 
shows the as—c^s marginalized distribution, for three data 
combinations. We note that as is significantly correlated 
with other cosmological parameters, in particular with 
ns(fco), Tc and as, so applying further priors to as moti- 
vated by inflation theory would affect these results, but 



it is useful to see what the data imply without such im- 
positions. 

4.3. Constraints on as from the High £ Excess Power 

Intrinsic CMB anisotropics on small angular scales are 
expected to be significantly suppressed by photon viscos- 
ity and the finite thickness of the last scattering region. 
Data from the CBI were the first to show this damping 
tail by mapping a drop of more than a factor of ten in 
power between £ = 400 and £ = 2000. This damping has 
subsequently also been observed by both ACBAR and 
the VSA. 

A number of effects are expected to produce sec- 
ondary anisotropies which peak at high £ . These 
include the Vishniac effect fVishni a^ 119871) . gravita- 
tional lensing (Blanchard & Schneiderl ll987t) . patchv re- 
ionization (Ae hanim et al.lll996fl . the Sunyaev-Zeldovich 
effect in galaxy clusters at moderate r edshifts z < 5 
fe.g.. lBond fc MvergllToM lCooravll200H) . and Sunyaev- 
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Fig. 8. — CMB constraints on Qi,h^ and ris, marginalized over 
other parameters. Shown are the la and 2(t constraints from: 
"WMAP only" (dashed lines delineating the blue region); "WMAP 
+ CBI" (dash triple-dot lines delineating the orange region); "CBI 
+ ALL" (dash-dot lines delineating the green regions) . In all cases 
a flat plus weak-h prior is used. 



0.50 




Fig. 9. — CMB constraints on tq and erg, marginalized over other 
parameters. Shown are the Icr and 2a constraints from: "WMAP 
only" (dashed lines delineating the blue region); "WMAP + CBI" 
(dash triple-dot lines delineating the orange region); "CBI -|- ALL" 
(dash-dot lines delineating the green regions). In all cases a flat 
plus weak-/i and LSS priors are included. Only data at £ < 2000 
is used in this analysis, although the trg results are consistent with 
constraints from the high-£ analysis of secondary SZE anisotropy. 



Zeldovich fluctuations from the first stars at higli red- 
sfiifts (z ~ 20) (Oh, Coorav. & Kamionkowski 2003). 

We previously considered the possible implications 
of the SZE in g alaxy clusters at moderate redshifts 
ijBond et alJl2004j) . and here we discuss this effect in the 
light of the new results presented above. We have es- 
timated cTg by fitting jointly for a primary CMB spec- 
trum an d template SZE spectr a. The technique is de- 
tailed in lGoldstein et a l] (I2003D where a combination of 
high-f bandpowers ijMason et alJl2003t iKuo et"al]l2004t 
IDawson et al.ll2002|) was used in a two parameter fit of 
"primary" and "secondary" spectrum amplitude param- 
eters. The SZE contribution to the power spectrum is 
highly dependent on the amplitu de of the mass fiuctu- 
ationS j characterised by an^ (e.^-- jKomatsu fc: Kitavamal 
IT999t iSeliak et aLllMl IBond et all 120041) . Since the 
SZE power spectrum has a weak dependence on flj, in 
addition to a strong cg dependence, it is useful to use 
an amplitude parameter crf^ to describe the scaling of 
the secondary SZE power spectrum. Assuming that 




Fig. 10. — Marginalized distributions for the running spectral 
index parameter Os defined in equation |5] The blue/dashed curve 
is for "WMAP only" and red/dotted is for "CBI + WMAP' (both 
for the flat plus weak-h prior case). The green/solid curve shows 
the result when also including our LSS prior. 
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Fig. 11. — Marginalized 2-d distributions for the running spectral 
index parameter and amplitude erg. Shown are the la and 
2a constraints from: "WMAP only" (dashed lines delineating the 
blue region); "WMAP + CBI" (dash triple-dot lines delineating 
the orange region) both for the flat plus weak-h prior case. The 
dash-dot lines delineating the green regions show the results for 
the "CBI -f ALL" case with the LSS prior also included. As in 
Figure|n] only data at ^ < 2000 are used. 



the power spectrum Cf^ scales as (r2h/i)^(T|, we define 
crf^ EE (f]{,/i/0. 035)2/'^ tJs. It should be noted that sec- 
ondary anisotropics, unlike intrinsic anisotropics, are not 
expected to have a Gaussian distribution. Although the 
detections in these bands are marginal, the strong depen- 
dence of the SZE power spectrum on the linear ampli- 
tude of the matter power spectrum already implies some 
constraints on the value of ag. The primary spectrum 
amplitude parameter encompasses the linear amplitude 
of perturbations as well as the effects of n., and tc on 
the expected high-^ bandpower. IGoldstein et al.l 1)20031 ) 
present an extensive discussion of the fitting procedure. 

The method approximates the effect of the non- 
Gaussian secondary anisotropy power spectra by mul- 
tiplying the expected sample variance in each band by 
a factor /ng of between 1 and 4. The bin covariances 
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are increased by the same factor. While this approach 
is simphstic, it is supported by numerical simulations 
which have shown the variance of simulated power spec- 
tra to be greater than the Gaussian case by a factor of 
approxim ately 3 for the £-rang e considered fCoorav 200 
[white Hernauist. fc SnringeJ 120021: iKomatsu fc Se li 
120021 iZhang. Pen, k Wand 12002(1 . Future work may 
require a more accurate treatment of non-Gaussianity. 
However we note that even large changes (/ng = 1~4) 
have a minimal impact on the results. There are also 
theoretical uncertainties of a factor of 2 in the theoret- 
ical SZE power spectrum predictions. These theoretical 
uncertainties translate into ^ 10% in erg and are also a 
limiting factor. 

For this work we used the last two bands of the 
power spectrum in the "even" binning of Table ^ for 
the C BI results, the la st three bands of the ACBAR re- 
sults l|Kuo et aLll2004D. and the two band result from 
the BIMA arrav l|Dawson et all 120021) . As a template 
primary spectrum we used the best-fit ACDM model 
with power law initial spectrum to the WMAP data^*. 
We assign a Gaussian prior with an rms of 10% around 
the best-fit amplitude for the primary spectrum while 
keeping all other parameters fixed, and marginalize over 
the primary amplitude parameter when deriving the final 
confidence intervals for cr|^. We have also included, for 
the CBI bandpowers, uncertainties due to the residual 
discrete source and thermal noise corrections. By con- 
sidering the of the CBI-fACBAR-HBIMA to a model 
comprising primary anisotropy and zero SZE signal, and 
with /ng = 1; we associate a statistical significance of 
98% with the detection of an SZE foreground at £ > 2000. 
The BIMA+CBI data alone give a 92% significance. 

In Table El we show the constraints on crf^ obtained 
from the fits to CBI -t- BIMA, and to CBI + BIMA -f 
ACBAR. The distributions have long tails extending to 
low values of cr|^ and are effectively unbounded from 
below (see Figure I12|l . In the context of our calcula- 
tion this is entirely due to Gaussian statistics and the 
results of changing variables to bandpower^^^ (in effect). 
We therefore define the confidence intervals as centered 
around the maximum in the distribution with the 1-a 
bounds given by a drop of a factor of e~^/^ on either 
side. 

We note that the results we derive for erf ^ are rather 
similar to those obtained using the one-year deep CBI 
fie ld in conjunction with B IMA and ACBAR, as reported 
bv lGoldstein et all l|2003() . We have repeated this analy- 
sis of the CBI one-year deep field results using the cross- 
calibration with WMAP, and find similar results. This is 
because the deep field component of our combined two- 
year data dominates the high £ power, and this is not 
changed much when the extra deep field 2001 data are 
added. What is important to note is that the much larger 
spatial coverage afforded by the full mosaic dataset (and 
thereby lesser sample variance) does not much diminish 
the amplitude of &. 
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5. CONCLUSIONS 



The CBI power spectrum is compared with WMAP 
and ACBAR results in FigurelTSI These results, together 
with those from a host of other ground- and balloon- 

^® |http : //lambda ■ gsf c ■ nasa ■ gov | 




Fig. 12. — Constraints on secondary SZE anisotropy, as param- 
eterized by the effective parameter erg . The curves are obtained 
from fits to the data at £ > 2000 (CBI and BIMA (red dashed), 
and CBI, BIMA, and ACBAR (blue soUd)). The fitting accounts 
for a separate contributions from template primary and secondary 
spectra. The marginalized distribution is heavily skewed towards 
low cr|^ values due to the assumed scaling of the secondary sig- 
nal. In the inset we have plotted the distribution against 
to show how the distribution is approximately Gaussian in this 
variable, which roughly corresponds to the high-£ bandpower in a 
given experiment. 



based experiments in recent years, are consistent with the 
key predictions of structure formation and inflationary 
theories: The universe is close to fiat; the initial spectrum 
of perturbations is nearly scale invariant; oscillations and 
damping in the power spectrum evince the expected sig- 
natures of sub-horizon scale causal processes; initial con- 
ditions are Gaussian, and are consistent with adiabatic 
fiuctuations; and the magnitude of fiuctuations from the 
largest scales down to galaxy cluster scales is consis- 
tent with what is needed to produce locally observed 
structures through gra vitatio nal col lapse. For discus- 
sion o f these points see : Bond et al.l l)2002D . iPeiris et"al] 
l)2003tl . and references therein. The concordance of ob- 
servational results with theoretical expectations has per- 
mitted cosmological parameters to be determined with 
precision. In this work we obtain: flbh"^ = 0.0225] 



n,h'' = O.Ultroll = 0-74^''oM, rc = O.ll^U 
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from a selection of current primary anisotropy data in- 
cluding CBI, WMAP, ACBAR, and Boomerang, and us- 
ing the fiat plus weak-h prior (see Table Similar re- 
sults are obtained when large-scale structure priors are 
incorporated (Table 2J) . As discussed in § 0] a flat prior 
(i.e., assumption that Otot = 1) is imposed on most 
of our parameter analysis; while supported by observa- 
tional data this does impose a strong constraint, and 
some parameter estimates would be less accurate with- 
out it. A marginal detection of a running scalar spectral 
index remains, and is consistent with that presented by 
iSpergel et al (20031. 

As discussed in § 14.11 the addition of CMB data from 
600 < i < 2000 significantly improves constraints on 
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Fig. 13. — The CBI+ M^MAP+ ACBAR Spectrum. The solid black line is the WMAP ACDM model with a pure power-law primordial 
spectrum (wmap_lcdm_pl_inod6l_yrl_vl . txt).The highest-^ ACBAR, point has been displaced slightly to lower £ for clarity. 



TABLE 6 
Values of (t|^ 







/ng = 1 


/ng = 2 


/ng = 3 


/ng =4 


CBI H 
CBI H 


- BIMA 

- BIMA H 


h ACBAR 0.981°;°^ 


0.98«:»? 


0.981H? 


0.96ini 
0.9810;°^ 



Note. — (t|^ values derived from the marginalized distributions obtained 
by fitting an SZE spectrum to the high-£ CMB data. The value for /ng is 
the factor used in rescaling the sample variance for each band (and inter- 
band correlations) to approximate varying degrees of non-Gaussianity. We 
find that the confidence limits do not depend strongly on the assumed /ng- 



i7f,/i^, Us, the amplitude of the prhiiary anisotropies, 
the age of the universe, and tc relative to what is ob- 
tained with only large-scale CMB data (see Figure d) . 
In the absence of a restrictive tc prior the £ < 600 
data leave significant degeneracies which are broken by 
the higher-^ experiments (see Figures |S1 and E))- We 
note that the improvement between the "CBI-I- WMAF^ 
and "CBI-I-ALL" cases comes primarily from the ad- 
dition of the Boomerang data. Improvements are also 
seen in analyses which allow a running scalar spec- 
tral index (Figure The tight constraint on the 
baryon density, Q^h'^ — 0.0225j;o;QOQg compares favor- 
ably with observationally determi ned BBN va lues of 
Qbh^ = 0.0214 ± 0.0020 ( Kirkmai TeTaLl l200l . We 
have also obtained an accurate measurement of Ug from 
the CMB data only, = 0.95 ± 0.02. These resuhs 
are robust with respect to prior assumptions, such as 
flatness, imposed on the analysis. By way of compar- 



ison the WMAP-only values for these parameters are 
Qbh^ = 0.0243 ± 0.0016 and = 1.01 ± 0.05. The 
breaking of these degeneracies largely relies on the ratio 
of power levels on small angular scales to those on large 
angular scales, so the precision of these results has ben- 
efited from the accurate cross-calibration with WMAP. 
The CBI data also favor a negative running scalar spec- 
tral index a, = -0.087 ± 0.028 (CBI-hALL-hLSS), con- 
sistent with the results from WMAP combined with LSS 
constraints 

In Figure [TH we show the same data as plotted in Fig- 
ure^] now on a log- log plot and with additional curves 
which show the expected level of SZE power for t he two 
sets of simulations discussed bv lBond et alJ l)2004j) . Note 
that the fortuitous "agreement" between the CBI and 
ACBAR power levels at £ > 2000 is not expected if the 
power has a significant component due to the Sunyaev- 
Zel'dovich Effect because of the different observing fre- 
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Fig. 14. — The CBI+ WMAP+ ACBAR Spectrum + high £ points from BIMA. The curves at high £ show the levels of SZ power expected 
in representative models using moving mesh hydrodynamics simulations (dotted) and smooth particle hydrodynamics (dashed) simulations 
(see text). The green and pink curves correspond to 30 GHz and 150 GHz, respectively. In these simulations (t|^ = 0.98, which also fits 
well the WMAP and CBI observations at lower £ for the case of a running spectral index (see Table 5). The highest-^ ACBAR point has 
been displaced slightly to lower £ for clarity. 



quencies. Nevertheless, given the uncertainties in these 
two measurements, it can be seen that the models span 
a range of power at high £ which fits both the CBI and 
ACBAR observations. 

The detection of power at ^ > 2000 i s cons istent with 
the results presented bv IMason et al.1 lj2003|) . although 
somewhat lower. We find a bandpower 355t}22MK^ 
(68% confidence, including systematic contributions). 
By combining this result with high-£ results from BIMA 
and ACBAR we detect power in excess of that expected 
from primary anisotropy at 98% confidence. This re- 
sult includes a marginalization over expected primary 
anisotropy power levels. Assuming the signal in excess 
of expected primary anisotropy is due to a secondary SZ 
foreground we determine crf^ = 0.96t^:EJ? (68%). The 
lower confidence level of the detection of an excess, and 
also the smaller values of crf^, are chiefly due to the 
lower high-€ bandpower we obtain and the inclusion of 
the uncertainty in the primary anisotropy bandpower 
at £ > 2000. The strong dependence of the observable 
power on as gives rise to firm upper limits on as but a 
tail to low values (Figure lT^ . It should be borne in mind 
that there are systematic uncertainties in the theoretical 
prediction of the power spectrum due to secondary SZ 



anisotropics which correspond to a 10% systematic un- 
certainty in as- 

An appreciable fraction of CBI data were rejected by 
vetoing NVSS sources, and furthermore the uncertainty 
in the power level of the source population remaining 
after the NVSS veto is a limiting factor at i > 2000. In 
late 2004 a sensitive, wideband continuum receiver will 
be commissioned on the Green Bank Telescope (GBT) 
to deal with both of these issues. This will result in a 
more sensitive determination of the total intensity power 
spectrum at all i covered by the CBI. Since the end of 
the observations reported here, the CBI was upgraded 
and dedicated to full-time polarization observations. 
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